VANISHING VISCOSITY AND THE ACCUMULATION OF 
VORTICITY ON THE BOUNDARY 



JAMES P. KELLIHER 

Abstract. We say that the vanishing viscosity hmit holds in the clas- 
sical sense if the velocity for a solution to the Navier-Stokes equations 
converges in the energy norm uniformly in time to the velocity for a 
solution to the Euler equations. We prove, for a bounded domain in 
dimension 2 or higher, that the vanishing viscosity limit holds in the 
classical sense if and only if a vortex sheet forms on the boundary. 



As is well known, for radially symmetric initial vorticity in a disk, the velocity 
for a solution to the Navier-Stokes equations converges in the energy norm 
uniformly in time to the velocity for a solution to the Euler equations. It 
was shown recently in [8j that for such initial data in a disk, it also happens 
that a vortex sheet forms on the boundary as the viscosity vanishes. By 
a vortex sheet, we mean a velocity field whose vorticity, as a finite Borel 
signed measure (an element of the dual space of C{^})), is supported along 
a curve — the boundary, in this case. 

It turns out that this phenomenon is in a sense more universal: either 
both types of limits hold or neither holds for an arbitrary bounded domain 
in R"', d > 2, with C^-boundary and with no particular assumption on the 
symmetry of the initial data. More precisely, the vanishing viscosity limit in 
the classical sense (condition (B) of Section [2]) holds if and only if a vortex 
sheet of a particular type forms on the boundary (conditions (E) and (-E2) 
of Section [2]). Now, however, the vortex sheet has vorticity belonging to the 
dual space of H^{Q) rather than C(r2). We show this in Theorem 12.11 for 
no-slip boundary conditions and in Theorem 14.11 for characteristic boundary 
conditions on the velocity. 

1. Introduction 

Let be a bounded domain in M'^, d > 2, with C^-boundary F, and let n be 
the outward unit normal vector to L. A classical solution {u,p) to the Euler 
equations satisfies 

, . ( diJl + u ■ Vn + Vp = f and divn = on [0, T] x Q, 
^ ^ \ u- n = on [0,r] X r andu = n° on {0} X f). 
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These equations describe the motion of an incompressible fluid of constant 
density and zero viscosity. The initial velocity TiP must at least lie in 

H = ^ue {L^{^)y : divti = in 0, u • n = on r| 

endowed with the L^-norm, which, along with 

y = |n G {H''{n)y : divu = in 0, -u = on pj 

endowed with the i?^-norm, are the classical spaces of fluid mechanics. 

We assume that uP is in C^+'(Q) n H, e > 0, where /c = 1 for two 
dimensions and k = 2 for 3 and higher dimensions, and that / is in C^([0, t] x 
il) for all t > 0. Then as shown in [7] (Theorem 1 and the remarks on p. 
508-509), there is some T > for which there exists a unique solution, 

win C^([0,T];C'=+^(S7)), (1.1) 

to [EE). In two dimensions, T can be arbitrarily large, though it is only 
known that some nonzero T exists in three and higher dimensions. 

The Navier-Stokes equations describe the motion of an incompressible 
fluid of constant density and positive viscosity u. A classical solution to the 
Navier-Stokes equations with no-slip boundary conditions can be defined in 
analogy to {EE) by 

dtu + u ■ Vu + Vp = vAu + / and div u = on [0, T] x 0, 
u = on [0, T] X r and n = n° on {0} x n, 

where is in H and / is in L^{[0,T];L'^{n)). We wih work, however, with 
weak solutions to the Navier-Stokes equations. (See, for instance. Chapter 
III of [lOj.) Such weak solutions lie in L°°([0, T]; i/) n L^{[0,T];V). 

In Section [2] we prove various equivalent forms of the vanishing viscosity 
limit for no-slip boundary conditions, including the formation of a vortex 
sheet on the boundary, and remark briefly on their derivation in Section [3j 
In Section H] we extend the results of Section 2 to characteristic boundary 
conditions. We discuss, in Section [5l our results in relation to those in [8] 
on vortex sheet formation for a disk. Finally, in Section [6l we include some 
technical lemmas employed in Section [2j 

2. Equivalent forms of the vanishing viscosity limit 

Let u he a. classical solution to (EE) in 0, and u = u'^ he a weak solution 
to (NS) in Q as in Section [U and assume that ^TiP in H and / — s- / in 
L^{[0,T];L'^{n)) as i/^O. 

Let 7n be the boundary trace operator for the normal component of a 
vector field (see Lemma [6TT]) . Let M{Q) be the space of finite Borel signed 
measures on Q — A4(0) is the dual space of C{Q). Let /i in A4{Q) be the 
measure supported on T for which /i|r corresponds to Lebesgue measure on 
r (arc length for d = 2, area for d = 3, etc.). Then fi is also a member of 
H\n)'. 



(NS) 
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We define the vorticity uj{u) to be the d x d antisymmetric matrix 

u;{u) = ^[Vu-{Vuf]. (2.1) 

When working specifically in two dimensions, we can alternately define the 
vorticity as the scalar curl of u: 

uj{u) = diu"^ - d2U^ . (2.2) 

Letting uj = u{u) and U = uj{u), we define the following conditions: 

(A) u weakly in H uniformly on [0, T], 

(A') u —*u weakly in {L?(yL))'^ uniformly on [0,T], 

{B) u^uin L^{[Q,T];H), 

(C) Vu^Vu- {-/n-,ufi) in {{H^{n))'^'"^y uniformly on [0,r], 

(D) Vn S/u in {H"^ (n))'^'"^ uniformly on [0, T], 

(E) w^cJ-i(7n(---^),¥^> in ((//I (17))'^^'^)' uniformly on [0,r]. 

In conditions (C), (D), and (E), the convergence is in the weak* topology 
of the given spaces. In (C) and (E), {{H^{Q))'^^'^y is the dual space of 
{H\n))'^'"^; in (D), {H-^n)^'^ is the dual space of {H^in)^'^. Thus, 
condition (C) means that 

(Vn(t),M) ^ (Vu(t),M) - j{M ■ n)-u{t) in L°°([0,r]) 

for any M in {H^{i~l))'^^'^, condition (D) means that 

(VM(t),M) ^ {Vu{t),M) in L°°([0,r]) 
for any M in (iifg (Q))'^^'^, and condition (E) means that 

(u;(t), M) ^ {uJ{t), M) - i ^((M - M^) • n) • u{t) in L-([0, T]) 

for any M in (H^n))'^'"^. 

In two dimensions, defining the vorticity as in Equation (j2.2p . we also 
define the following two conditions: 

(£^2) uJ^iJ-{u-T)fim {H^in))' uniformly on [0, T], 

(-F2) a; — > uJ in H^^{0,) uniformly on [0,T]. 

Here, r is the unit tangent vector on T that is obtained by rotating the 
outward unit normal vector n counterclockwise by 90 degrees. 
Condition {E2) means that 



jm-T)fin L-([0,r]) 



(u;(t),/)^(uJ(t),/) 

for any / in H^(Q), while condition {F2) means that 
(u;(t),/)^(cJ(t),/) in L^{[0,T]) 



4 JAMES P. KELLIHER 

for any / in Hq{Q). 

Theorem 2.1. Conditions (A), (A'), (B), (C), (D), and (E) are equiv- 
alent. In two dimensions, conditions (E2) and (F2) are equivalent to the 
other conditions. 

Proof. (A) (A'): Let v be in By Lemma [Ol v = w + Vp, 

where w is in H and p is in H^{Q). Then assuming (A) holds, 

{u{t),v) = {uit),w) ^ {u{t),w) = {uit),v) 

uniformly over t in [0,T], so (A') holds. The converse is immediate. 

(A) <^=^> (B): The forward implication is proved in Theorem 1 of [3]. The 
backward implication is immediate. 

(A') =^ (C): Assume that {A') holds and let M be in (H^n))'^'''^. Then 
{Vu{t),M) = -(n(t),divM) ^ -(u(t),divM) in L°°([0,r]). 

But, 

-(M(t),divM) = {Vu{t),M)- j (M-n) -u, 

giving (C). 

(C) =^ (D): This follows simply because H^{n) C H\n). 

(D) =^ (A): Assume (D) holds, and let v be in H. Then v = divM for 
some M in {H^{n))'^'"^ by Corollary [631 so 

{uit),v) = {u{t),divM) = -{Vu{t),M) + j {M-n)-u{t) 

= -{Vu{t),M) -{Vu{t),M), 
uniformly over [0, T]. But, 

-{Vu{t),M) = (n(t),divM) - j {M ■n)-u{t) = 

from which (A) follows. 

Now assume that d = 2. 
(A') =^ (E2): Assume that {A') holds and let / be in H^{Vt). Then 
{u{t)J) = -(divn^(t),/) = {u^{t),Vf) = -{u{t),V^f) 
^-(n(t),V^/) in L-([0,T]) 
where «-*- = — (u^, u^) and we used the identity u!{u) = — divn-*-. But, 

-{u{t),V^f) = (u^it), V/) = -(divn^(t), /) + l^iuit)^ . n)f 

= -(divll^(t),/) - ■ T)f = (uJ(t),/) - • T)f, 

giving {E2). 



VANISHING VISCOSITY AND VORTICITY ON THE BOUNDARY 5 

(E2) =^ (F2): Follows for the same reason that (C) =^ (D). 

(F2) =^ (A): Assume (-F2) holds, and let v be in H. Then v = V"*-/ for 
some / in i?g(Q) (/ is called the stream function for v), and 

{uit),v) = (n(t), V^/) = -(u^{t),Vf) = (divn^(t),/) 
= -(c^(t),/) - -{Ht)J) in L-([0,r]). 

But, 

-(clJ(t),/) = (divu^(t),/) = V/) = V^/) 

which shows that (A) holds. 

What we have shown so far is that (A), (A'), {B), (C), and (D) are 
equivalent, as are (-B2) and (F2) in two dimensions. It remains to show that 
{E) is equivalent to these conditions as well. We do this by establishing the 
implications (C) =^ (E) =^ (A). 

(C) =^ (E): Follows directly from Equation ()2.ip . 

(E) =^ (A): Let w be in and let X be the vector field in (i?2(j7)ni7^(0))'^ 
solving Ax = V on ^} (x exists and is unique by standard elliptic theory). 
Then, utilizing Lemma 16.61 twice (and suppressing the explicit dependence 
of u and u on t), 

(u, v) = (n. Ax) = -(Vti, Vx) + j (Vx • n) • u = -{Vu, Vx) 
= —2{Li!{u),uj{x)) — J (Vux) ■ n = —2{lj{u),uj{x)) 
-2{uj{u),Lo{x)) +2-1 ((lj(x) - uj{x f) • n) • u 



2 

-2(a;(?l), a;(x)) + 2 {uj{x) -n) -u 



(2.3) 

= -(VH, Vx) + J (Vux) ■n + 2 j (u;(x) • n) • H 

= -(Vu, Vx) + 2 j {uj{x) • n) • u 

= (u. Ax) — j (Vx • n) • u + 2 y (<j-'(x) • n) • u 

= (n, v) — j ((Vx)"^ • n) • n. 
Thus, (£;) ^ (A) if and only if 

j ((Vx)^-n)-u = 0. (2.4) 
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But, (div(Vx)"^)-' = djdiX^ = didwx or div(V3;)-^ = Vdivx. Similarly, 
div(VTi)"^ = VdivM = 0. It follows that 

j ((Vx)^ • n) • u = ((Vx)^, Vu) + (Vdivx,tt) 

= (Vx, (Vu)^) — (divx, divn) + J u • ndivx = ((Vu)^, Vx) 
= -(div(Vu)^, x)+ J {{Vuf ■ n)x = 0. 

□ 

3. Remarks 

The equivalent conditions of Theorem 12.11 complement those of [1] , [11] , 
[H], and [6]. 

It is only in the proof of {A) =^ (B) — in which we quote a result of 
Kato's in [3] — where the requirement that n be a classical solution to the 
Euler equations and that f f in L^([0,T];L^) is used; in fact, it is the 
only place where the fact that u and u are solutions to the Navier-Stokes and 
Euler equations, respectively, appear in the proof at all. That is, assuming 
only that u is a vector field parameterized by i' that lies in L^{[0,T]; H) n 
L'^{[0,T];V) and that u is a vector field lying in L°°([0,T];F n H^{n)), 
all of the implications in the proof of Theorem 12.11 remain valid except for 
{A) =^ {B). 

The proof of (A) =^ (B) in [2] consists, using our terminology, of 
proving {B) =^ (A) =^ (i) =^ (ii) =^ (B), where (i) and {ii) are 
the conditions, 

rT 

(i) V \ ||Vn||^2m) ^ 0, 







T 



{it) V / ||Vn||^2(re,) ^ 0' 

J 

with Yf-v a boundary layer of width proportional to v. 

If we weaken the regularity of F from C"^ to only locally Lipschitz, then 
the proof of {E) =^ {A) fails because we would only have x in (//Q(il))'^. 
Kato's proof that {A) =^ (B) also requires a boundary. 

In two dimensions, we need only have convergence of the vorticity away 
from the boundary — condition (F2) — to insure that the vanishing viscosity 
limit holds. In particular, it follows that formation of a vortex sheet on the 
boundary of a type other than that given in (-E'2) is inconsistent with u being 
a solution to (NS). In higher dimensions it is an open problem whether the 
analogous statement is true; that is, whether (F) =^ (A), where (F) is 
the condition, 

(F) LO^uJ in {H-'^in)^'^ uniformly on [0,r]. 
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The remarks that fohow attempt to give some insight into the nature of this 
problem. 

One approach to proving that (F) =^ (yl) is to prove that (F) =^ (D), 
since we have (D) =^ (A). So suppose that (F) holds, and let M be in 
{Hq{Q))'^^'^. For any vector field v, 

{Vv, M) = {Vv - (Vvf, M) + ((Vv)^, M) = 2{uj{v), M) + (Vv, M^). 
Thus, 

(V-u, M - M^) = 2{lo{u),M) 2{uj{u),M) = {Vu, M - M^). 

If M is antisymmetric then M — M'^ = 2M and we conclude that (D) holds 
for antisymmetric matrix fields in {Hq{Q))'^^'^ . But if M in (i^Q (O))'^^'^ is 
symmetric, 

2{uj{u),M) = {Vu, M) - {{Vuf, M) = {Vu, M - M'^) = 0, 

so {oj{u),M) = {uj{u),M) = 0, and we can conclude nothing from this 
approach. 

But some use can still be made of this observation. Let v be any element 
of H. Then from Corollarv 16.51 for some M in {Hq{^}))'^^'^, 

{u,v) = (u,divM) = -(Vn,M). 

Now, if we could insure that M can be chosen to be antisymmetric, then if 
(F) holds for M so does {D), as we just showed, and 

-(Vn,M) -{Vu,M) = {u,v), 

and {A) would follow. 

In two dimensions, we can choose 

M=(° "/), (3.1) 

where / is the stream function for v as in the proof of {F2) =^ {A), which 
gives a slight variation on the proof of that same implication. 

In three dimensions, however, it is not possible to find such an M. To see 
this, suppose that M in {Hq{Q))'^^'^ is antisymmetric. Then can write M as 

/ a b 
M = i-a c 
\-b -c 

with a = 6 = c = 0onr. In this form the condition divii = divdivM = 
is automatically satisfied, and letting cD be the vector (c, —b, a), we see that 

li = divM = curia), (3.2) 

where curl is the usual three-dimensional operator. But curl maps Hr]C°^{^}) 
bijectively onto itself when T is (see, for instance, [2]), so in general we 
only have cu • n = on F. That is, the condition that M be antisymmetric 
is not compatible with the condition that it vanish on F. 
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Finally, let 

E{n) = |u G divn G L'^{n),u- n = on r| , 

with ||'u||£;(f7) = ||'^|lL2(f^) + ||divti|[j;^2(Q). It is easy to see from the proofs 
of {A') (C) and (D) (A) that condition (D) can be weakened 

from convergence in the dual space of (//q (17))'^^'^ to convergence in the dual 
space of {E{Q))'^. (This is advantageous as a sufficient condition, though not 
as a necessary one.) 

Returning to Equation (13. 2p . the condition ui n = on F does not translate 
to M ■ n = on r. Hence, M does not lie in {E(Q,))'^ so we cannot use this 
weakening of condition (D) to conclude that {A) holds. 

4. Characteristic boundary conditions on the velocity 

We modify (NS) by allowing the velocity on the boundary to be equal 
to a nonzero time-varying vector field 6, which is required, however, to be 
tangential to the boundary. This gives 

(ArQ\ / ~^ ^ ' ~'~ ~ + / ^^'^ '^^^ u = on [0, T] x Q,, 
l^V^bj I n = 6 on [0,r] X r and u = < on {0} x Q, 

where, as before, is in H. 

We require sufficient regularity on b so that (NSb) is well-posed. For 
simplicity, we will assume that 6 • n = on F with 

6g L~([0,T];/73/2(r)), dtbe L^{[0,T];H-^/\T)) (4.1) 

so that b lifts (extends) to a vector field, which we also call b, with 

b G L°°([0, T];Hn H^in)), dtb G L^qq, T]; ^^(J^)). 

The assumption on b in Equation (j4.ip is not the weakest possible, but the 
assumption on dfb can probably not be weakened. 

We can then use the equation that corresponds to n — 6 (which lies in 
L°°([0, T];H)r) L^([0, T];V) for classical solutions) to define a weak solution 
to (NSb). (This is essentially what is done in Section 4 of |12j.) 

With such solutions to (NSb) in place of those for (NS) — but without 
changing the formulation of {EE) — we define the condition, 

{C'') Vu^Vu- {-fn-,{u-b)fi) in ((//^(l^))'^^'^)' uniformly on [0,r], 

^E') ^^uJ-l (^„(. _ .^), (u - 6)^> in i{H\n)r>'''y 

uniformly on [0,T], 
and in two dimensions, the condition, 

(^2) co^uj-{{u-b)-T)fi in {H^in))' uniformly on [0,T]. 
Theorem 12.11 then becomes: 
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Theorem 4.1. Let u he a solution to (NSh) andu be a solution to (EE). 
Conditions (A), (A'), (B), (C^), (D), and (E^) are equivalent. In two 
dimensions, conditions (E2) and (F2) are equivalent to the other conditions. 

Proof. The proof of [A') =^ (C^) is identical to the proof of {A') =^ (C) 
except that a boundary term is included in the first step: this term leads 
to the "—6" in condition (C*). A similar comment applies to the proof of 
(A') =^ {El) and (C7) =^ {E^). The proof of (C) =^ {D) and 
(£'2) =^ {E2) are unaffected by the presence of the vector field 6, 
all other implications except for [A) =^ [B) and {E^) =^ {A). 

(A) =^ (B): In [12] it is shown, using our terminology, that {B) =^ 
(i) =^ (ii') =^ (5), where {ii') is the condition 

rT rT 

{ii') V ||Vrnr||i2(r,^^^j) ^Qoiu || Vrnn||i2(r^^^^) ^ 0. 

Here, V-7- is the gradient only in the direction tangential to the boundary, 
ut and Un are the components of the velocity tangential and normal to the 
boundary, respectively, and is a boundary layer whose width 5{v) is of 
arbitrary order larger than v. 

But, {B) =^ {A) is immediate, and {A) =^ {i) follows from com- 
bining the argument on pages 232 through 233 of [12j with the proof of the 
implication {A) =^ {i) on page 90 of ^ (in Kato's terminology, this is (ii) 
implies (iii)). This gives {A) =^ {B). 

(E^) =^ (A): Adapting the argument of {E) =^ {A), we see that the 
first boundary integral in Equation (j2.3p does not vanish, since now u = b 
on r. Also, u becomes u — b in the boundary integrals involving uj{x). This 
leads to 

{u, v) = {u, Ax) + j (Vx -n) - b- j (Vx • n) • u + 2 j {uj{x) ■ n) ■ {u - b) 

= {u, v)- j ((Vx)^ • n) • (n - b). 

The last boundary integral vanishes for the same reason that Equation (j2.4p 
holds, u — b being in H, completing the proof. □ 

5. Radially symmetric initial vorticity in a disk 

We assume, in this section only, that is the unit disk D and that the 
initial vorticity is radially symmetric. In this case, the solution to {EE) 
is stationary: uj{t) = uf' for all time t. 

The vanishing viscosity limit in the classical sense of condition {B) holds 
in this setting under fairly general circumstances. Under the assumptions of 
Section [T] on the regularity of the initial velocity, and assuming that 6 = 0, 
this convergence is implicit in [4] at least for zero forcing (see [5]), but was 
first explicitly proved by Matsui (see [9j). For nonzero b having the regu- 
larity assumed in Equation ()4.ip . the convergence is a simple consequence 
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of the condition (ii') in the proof of Theorem 14.11 as observed by Wang in 
|12j . For substantially lower regularity on uP and on b than we assume, the 
convergence is established in [8]. 

More precisely, the authors of \8] assume that = vP and / = / = 
(which are not significant limitations, since one can handle u*^ — > ti*^ in 
H by using the triangle inequality, and nonzero forcing presents no real 
difficulties), with 

n° G R^{D) = jt; G {L^{D)f: v{x) = s{\x\)x^ for some s, uj{v) G L^-D)} 

= {f G H : u){v) G L^{D), uj{v) radially symmetric} . 

They assume that b{t, •) = a{t) — that is, b{t) is constant on the boundary — 
and that a G BV([0,r]), the space of bounded variation functions. They 
prove (combining Propositions 9.6 and 9.7 of [8]) that 

uj ^ uj - {B {2TTy^ -b-T)fi in M(D) uniformly on [0,T], (5.1) 

where 

B = / aJ°. 
Jd 

But, on T, ■ T is constant, so by Green's theorem, 

B = jifi ■T = 27ru°(x) • t{x) 

for any point x on T, and we see that Equation ()5.ip is the same as condition 
(£"2), except that the convergence is stronger. 

That is, both conditions (B) and (-©2) hold for a disk, except that the 
convergence in {E2) is in which is stronger convergence than that of 

(£'2)- What we have shown is that either both conditions {B) and {E2) hold 
or neither condition holds for a given initial velocity in a general bounded 
domain in the plane — and in the analogous sense, in M"^. It was the question 
of whether this was, in fact, the case that motivated this paper. 

The regularity we assume in Equation (j4.ip corresponds to a lying in 
H^{[0,T]), which is considerably stronger than the assumption in [8j that a 
lie in BV([0,T]). And their assumption on the regularity of vP is far lower 
than our assumption that vP lies in C^~^''{Q) Ci H. Without the assump- 
tion of radial symmetry, however, it seems unlikely that one can weaken our 
assumptions in Equation (|4.ip on b to any significant degree, since these 
assumptions go to the heart of establishing the existence of the correspond- 
ing weak solutions of (NSb). Weakening the regularity assumptions on 1^ 
would seem equally impossible, since the boundedness of Vn on [0, T] x S7 is 
indispensable in Kato's argument showing that {A) =^ (i?). 

6. Some technical lemmas 

In this section we assume only that 17 is bounded and that F is locally 
Lipschitzian, which of course includes the case that T is C^. 
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The various integrations by parts that we make are justified by Lemma [6.H 
which is Theorem 1.2 p. 7 of [10] for locally Lipschitz domains. (Temam 
states the theorem for boundaries but the proof for locahy Lipschitz 
boundaries is the same, using a trace operator for Lipschitz boundaries in 
place of that for boundaries: see p. 117-119 of [3j, in particular, Theorem 
2.1 p. 119.) 

Lemma 6.1. Let 

with \\v\\^(^^-^ = \\v\\i2^Q-j + ||divu||^2(Q). There exists an extension of the 

trace operator jn- (C'o°(^))'^ ~^ C°°(r) defined by u i—f u ■ n on F to a 
continuous linear operator from E{^}) onto H^^/'^{T). The kernel of is 
the space Eq{Q) — the completion of C^{0,) in the E(^}) norm. For all u in 
E{n) and f in H^Q), 



(n,V/) + (div7x,/) = ^(n-n)/. (6.1) 



Lemma 6.2. Assume that u is in {V'{n))'^ with {u,v) = for all V in V. 
Then u = Vp for some p in P'(Q). If u is in (L^(J7))'^ then p is in H^{Q); 
if u is in H then p is in H^{il,) and Ap = 0. 

Proof For u in {V'{n) f see Proposition 1.1 p. 10 of [lO]. For u in {L'^{Q)y 
the result follows from a combination of Theorem 1.1 p. 107 and Remark 
4.1 p. 55 of [3] (also see Remark 1.4 p. 11 of [lO]). □ 

Lemma 6.3. For any u in (L^(i7))'^ there exists a unique v in H and p in 
H^(Q) such that u = v + Vp. 

Proof. This follows, for instance, from Theorem l.lp. 107of[3j, which holds 
for an arbitrary domain, along with Lemma |6.2[ □ 

Lemma 6.4. For any f in L'^{Q) and a in {H^'^{T)Y satisfying the com- 
patibility condition, 



[ f= [ an 

Jn Jr 



there exists a (non-unique) solution v in {H^{Q))'^ to divv = f in Q., v = a 
on V . 

Proof. This follows from Lemma 3.2 p. 126-127, Remark 3.3 p. 128-129, 
and Exercise 3.4 p. 131 of [3j (and see the comment on p. 67 of [T]). □ 

Corollary 6.5. For any v in H there exists a matrix-valued function M in 
{Hl{n)Y'"^ such that v = divM. 

Proof. Let vheinH and observe that 

v"^ = / V ■ Vxi = — diYVXi + I {v ■ n)xi = 0. 
n Jo. Jn Jr 
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Thus, we can apply Lemma 16.41 to each component v"^ using a = to obtain 
a vector in Hq{Q.) satisfying divtD* = on fi, tn* = on T. Forming a 
matrix- valued function whose rows are ^nP' , . . . ^w"^ gives M. □ 

Lemma 6.6. Assume that u is in {H'^{Q,))'^ or {C^{Q))^ with divu = and 
V is in {H^{^})Y. Then 

(Vu, Vv) = 2{uj{u),uj{v)) + / (Vuv) ■ n. 



Proof. Directly from Equation (j2.1|) . 

2a;(u) • u;iv) = ^(Vn - (Vuf) ■ (Vv - (Vvf) 

2 (6.2) 

= Vu-Vv - (Vn)^ • Vv. 

Since divu = 0, we have (Vn)-^ • Vv = djv'^diu^ = dj{v^diU^) = div(Viiw), 
so if u and v are both in {C°°{Q))'^ with divn = then 

uj{u) ■ uj{v) = / Vu-Vv - (Vu)'^ - Vv 



= / Vu-Vv - / div(Vn?;) = j Vu-Vv - j {Vuv) - n. 
Jo, Jn Jn Jr 

The result then follows by the density of C°°(ri) in H'^in), H'^{Q), and 

c\n). □ 
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